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ABSTRACT 

It has been shown based on radio variability arguments that some BALQSOs 
(broad absorption line quasars) are viewed along the polar axis (orthogonal to ac- 
cretion disk) in the recent article of Zhou et al. These arguments are based on the 
brightness temperature, TJ, exceeding 10 12 ° K which leads to the well-known in- 
verse Compton catastrophe unless the radio jet is relativistic and is viewed along 
its axis. In this letter, we expand the Zhou et al sample of polar BALQSOs to 
the entire SDSS DR5. In the process, we clarify a mistake in their calculation 
of brightness temperature. The expanded sample of high TJ, BALQSOS, has an 
inordinately large fraction of LoBALQSOs (low ionization BALQSOs). We con- 
sider this an important clue to understanding the nature of the polar BALQSOs. 
This is expected in the polar BALQSO analytical/numerical models of Punsly in 
which LoBALQSOs occur when the line of sight is very close to the polar axis, 
where the outflow density is the highest. 

Subject headings: (galaxies:) quasars: absorption lines — galaxies: jets — (galax- 
ies:) quasars: general — accretion, accretion disks — black hole physics 



1. Introduction 

About 15% - 20% of quasars show broad UV absorption lines (loosely defined as ab- 
sorbing gas that is blue shifted at least 5,000 km/s relative to the QSO rest frame and 
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displ a ying a spread in velo city of at least 2,000 km/s) (jWeymann 1 119971 : iHewett and Foltz 
20031 ; iReichard et al.ll2003l ). The implication is that large volumes of gas are expelled as 
a byproduct of accretion onto the the central supermassive black hole. Understanding the 
details of this ejection mechanism is a crucial step towards understanding the physics of the 
central engine of quasars. Although evolutionary processes might be related to BAL out- 
flows, it is widely believed that all radio quiet quasars have BAL flows, but the designation 
of BALQSO depends on whether the line of sight intersects the solid angle 

subtended by the outflow. The standard model of quas ars is one of a hot accretion flow onto 
a black hole and a surrounding torus of molecular gas (lAntonuccil Il993l ) . The BAL outflow 
can be an equatorial win d driven from t 
above the molecular gas, iMurrav et all ( 
gions of the accretion flow (jPunsly 



re lum inous disk that is viewed at low latitudes, just 



1999a 



1995), or a bipolar flow launched from the inner re- 



b|). The most general feature of 3-D simulations of 



accretion flows onto black holes is huge ejections of gas from the central vortex of the accre- 
tion fl ow along the polar axis which rivals the accreted mass flux for rapidly spinning black 
holes (IDe Villiers et al 1120051 ; lHawley and Krolik II2006I ). Thus, it is of profound theoretical 
interest to look for evidence of these endemic polar ejecta. 

BALQSOs are so distant that direct imaging of the BAL region is beyond the resolution 
of current optical telescopes. Thus, much of the discussion of BAL geometry is based on 
de ductive rea s oning . A novel idea for finding the orientation of BALQSOs was developed 



m l Zhou et all (120061 ). They used radio variability information to bound the size of the the 
radio emitting gas then deduced that the radio emission must be viewed close to the polar 
axis and emanate from a relativistic jet, thereby avoiding the well known inverse Compton 
catastrophe. They found BALQSOs that satisfy these conditions in SDSS DR3. We expand 
their methods to SDSS DR5 and apply these results to the physics of the BAL wind launching 
mechanism. 



2. The Brightness Temperature 



Equation (5) of IZhou et all (120061 ) is not well defined since it is unclear what frame of 
reference is used to evaluate the quantities contained within. Apparently, this gave rise to 
estimates that are (1 + z) s larger than what we find in equation (2.4), below. In this section, 
we derive a formula for Physically, it is 7& evaluated in the cosmological frame of reference 
of the quasar, (T b ) q , that is relevant for assessing the "inverse Compton catastrophe." We 
want to express this in terms of observable quantities at earth designated by the subscript 
"o." First of all, the brightness temperature is the equivalent blackbody temperature of the 
radiation assuming one is in the Planck regime, hv <C fc&T. Consider a source in which the 
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monochromatic intensity has increased by an amount I(v) q in a time At q . The brightness 
temperature associated with the change in luminosity is 



bjg 



m^ 2 

1k h v\ 



From the monochromatic version of Liouville's theorem, iGunnl (119781 ) 



{FX 



= I(u) 



(2-1) 



(2-2) 



{F u ) , is the flux density observed at earth. The solid angle subtended by the source, Q Q , is 
bounded by the causality requirement that the source co uld not have expanded more rapidly 
than the speed of light during a time, At g , IGunnl ( 119781 ). 



a. 



proper area perpendicular to line of sight (c7rAt 



< 



(2-3) 



(angular diameter distance) 2 

where the angular diameter distance is (La- In a a cosmology with H$= 70 km/s/Mpc, Q\ = 
0.7 and Q m = 0.3, we can use the expression for cIa given by iPenl (119991 ). that is accurate to 
< 1% relative error, along with (2.1) -(2.3) to find: 



.0 x 10 12 (1 



b)q 



^ /lGHz) 2 (At /lyr) : 



r Z 2 (AF,(mJy)) °K 



(2-4a) 



Z = 3.31 - (3.65) 



x ([(1 + zf - 0.203(1 + zf + 0.749(1 + zf + 0.444(1 + z) + 0.205] °' 125 ) (2-4b) 

where Ai^(mJy) is the change in flux density in mJy measured at earth at frequency v Q 
during the time interval At a . 



When (Tb) q > 10 ° K, the inverse Compton catastrophe occurs. Most of the electron 
energy is radiated in the inverse Compton regime. Th e radio synchrotron spectrum from 
the jet is diminished in intensity to unobservable levels (jKellermann &: Pauliny-Tothl 119691 ). 
In order to explain the observed radio synchrotron jet in such sources, Doppler boosting 
is customarily invoked to resolve the paradox. Recall that for an unresolved source the 
observed flux density, (F v ) , is Doppler enhanced relative to the intrinsic flux density, (F v ) = 
(^"(F^intrinsic, where the Doppler factor, 5 is given in terms of T, the Lorentz factor of the 
outflow; 0, the three velocit y of the outflow; and a, is the spectral index, 5 = 1/[T(1— (3 cos^)] 
(ILind and Blandfordlll9851 ). Thus, Doppler beaming can cause (T b ) q to be enhanced by a 
factor of 5 3+a in (2.4). Assuming a flat spectral index for the unresolved relativistic jet, 
we choose a = 0. This allows us to find a minimum Doppler factor that avoids the inverse 
Compton catastrophe, S m i n [(Tb) q ] = [{Tb) q /(10 12 ° K)] 0,333 . If the jet plasma propagates with 
8 > fimin[{Tb)q)} then the intrinsic T& in the frame of reference of the jet plasma is sufficiently 
low that there is no inverse Compton catastrophe. 
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Table 1: Properties of the Radio Variable BAL Quasars a 



Name 


z Sp-.d&te Sedate 


17 


BI CIV 


bi Aim 


BI Mgll 


Type 


(SDSS J) 




(mJy) 


(mJy) 




(km/s) 


(km/s) 


(km/s) 




081102. 91+500724.5 t 


1.84 


23.07: 05/23/97 


18.80: 11/15/93 


5.7 


1568 


468 




LoBAL 


081618.99+482328.4 


3.57 


72.29: 05/01/97 


61.10: 11/15/93 


5.2 


2576 






HiBAL 


081839.00+313100.1 


2.37 


8.82: 10/23/95 


7.2: 12/15/93 


3.1 


6445 






HiBAL 


082817.25+371853. 7 C 


1.35 


21.18: 07/23/94 


14.50: 12/15/93 


10.5 




428 


1564 


LoBAL 


093348.37+313335.2 


2.60 


18.35: 10/23/95 


15.90: 12/15/93 


3.7 


4107 


1405 




LoBAL 


104106.05+144417.4 


3.01 


27.46: 12/99 


19.00: 12/06/93 


11.4 


2527 






HiBAL 


113445.83+431858.0 


2.18 


27.38: 2/20/97 


24.90: 11/15/93 


3.4 


10643 




1149 


LoBAL 


134652.72+392411.8 


2.47 


3.60: 08/19/94 


2.20: 04/16/95 


3.0 


1505 


557 




LoBAL 


142610.59+441124.0 


2.68 


6.78: 03/27/97 


5.0: 03/12/95 


3.6 


8652 


1389 




LoBAL 


145926.33+493136.8 d 


2.37 


5.22: 04/17/97 


3.60: 03/12/95 


3.4 


9039 


329 




LoBAL 


155633.77+351757.3 


1.49 


30.92: 07/03/94 


26.90: 04/16/95 


4.3 


e 






FeLoBAL 


165543.24+394519.9 


1.75 


10.15: 08/19/94 


8.50: 04/16/95 


3.0 


5805 






HiBAL 



a The first column in table 1 is the source name, followed by the redshift. Columns (3) and (4) are the FIRST 
and NVSS peak flux densities at 1.4 GHz, respectively. Column (5) is the statist ical significance of the r adio 
variability computed using (3.1). Columns (6) - (8) are the BALnicity indices. IWevmann et al.l (|199ll ). for 
CIV, A1II and Mgll, respectively. The last column denotes whether the source is a low ionization or a high 
ionization BALQSO. 

b Thcre is a new spectrum in SDSS with much higher S/N than the spectrum in lZhou et all (|2006l ) that clearly 
shows broad Al III absorption 

c from Zho u et a however there is now a much higher signal to noise spectrum in SDSS than the one 

used in lZhou et all ( 2006 ) and we have recomputed the BALnicity indices 
d from lZhou et all (|2006h 



identified as FeLoBALQSO in lBecker et all (|l997l ) 
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3. Polar BALQSOs 



Table 1 is a list of BALQSOs in which the radio variability requires that the jet is 
propagating well within 35 ° to the line of sight in order for the jet plasma to satisfy (TJ>) 9 < 
10 12 ° K. Since rather small variations of flux density create these conditions we must first 
prove that the sources are truly variable. We choose the condition for variability to be 

Sfp — Snp 



yf<j 



FP 



> +3 , 



(3-1; 



NP 



where Sfp and Snp denote the peak flux density at 1.4 GHz measured by the FIRST and 
NVSS surveys, respectively, <jfp and a^sp are the FIRST and NVSS peak flux uncertainties, 
respectively. The beam for the NVSS survey is considerably larger than FIRST (FWHM ~ 
45" versus FWHM ~ 5"). Thus, a compact source <C 5" (like a compact radio core which is 
the putative site of variable activity) should be detected with equal sensitivity if it exceeds 
the flux limit of the samples. The NVSS beam should pick up more extended flux than 
FIRST. Variable BALQSO fields in NVSS with confusion from nearby sources were omitted 
from table 1. Thus, a FIRST peak flux density larger than an NVSS peak flux density by 3a 
(hence the plus sign on the RHS of eqn. (3.1) ) is true va r iabilit y and not an artifact 
of the different beam sizes. One improvement from IZhou et all (120061 ) is that we used the 
peak NVSS flux density as opposed to the integrated NVSS flux density. This is preferred for 
the comparison described in (3.1), since the variable radio flux is most likely coming from 
sub-arcsecond regions of the jet and this choice is less sensitive to low surface brightness 
noise that can be picked up in the integrated flux (ICondonlll9971 ). 



It is crucial to assess the significance of the 3a results in table 1. In order to clarify 
the statistics, we must delineate the super sample from which the sources in table 1 are 
drawn. The potential variable BALQSOs in DR5 that can be detected by equation (3.1) 
require two cond i tions. First of all, the sources must have a BALnicity index > as defined in 
Weymann et al.l (ll99ll ) for at least one line either, CIV, A1III, or Mgll. Secondly, Sfp > 3.55 
mJy The NVSS survey has a peak flux limit of 2.5 mJy in general, b ut occasionally as low 
as 2.2 mJy and a minimum measurement uncertainty of ~ 0.45 mJy ICondon et al I (119981 ). 
Thus, Sfp < 3.55 mJy could never satisfy eqn. (3.1). There are 116 BALQSOs that satisfy 
these criteria and they form the master sample. Eqn. (3 .1) is a one sided pro bability, so if 
the errors are Gaussian distributed as in I Condon! (119971 ); I Condon et al I (119981 ) then there is 
less than a 20% chance of expecting a variable BALQSO from this sample of 116. However, 
we found 20 variable BALQSOs (12 of which made the table because Tg was sufficiently 
high to restrict the line of sight to < 35°). Thus, based on Gaussian statistics our findings 
are very significant. However, if the errors in the tail of the distribution ar e not Gaussian , 
but arise from unforseen sources of measurement error, then the analysis of ICondonl (119971 ); 
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Condon et al I (119981 ) does not apply and we have no means to assess the significance of our 
results. Note that there are 4 sources in table 1 with a variability that is significant above 
the 5<T level. 

Our BALQSO identifications and the BALnicity indices in table 1 were derived from 
spectra in the SDSS DR5 public archive. The SDSS spectra of radio-variable BALs were 
retrieved from the SDSS database and were analyzed using the IRA^A softwa r e. Fir st, the 
spectrum was de-reddened using the Galactic extinction curve, ISchlegel et al.l (119981 ). then 
the wavelength scale was transformed from the observed to the source frame. The spectra 
were fitted in XSPEC with a power law plus mult iple gaussian model, including fits for both 
emission lines and absorption lines (|Arnaudlll996l ). All the model parameters were kept free. 
The best fit to the SDSS data was determined using \ 2 minimization. If an emission bump 
around 2500 A is present then the continuum fit in the CIV (1550 A) region was extrapolated 
to longer wavelength. The BALnicity indices quoted in table 1 are a consequence of the 
method of spectral fitting described above and other methods might produce different results. 
However, the exact BALnicity index is not critical to this discussion, the sources in table 1 
clearly show BALs and that is the essential point of relevance here. 

In table 2, we calculate (T&) g , using equation (2.4) and the peak fluxes from table 1. 
From this, we calculate S m i n [(Tb) q ] in column (3) as described at the end of section 2. For 
each value of 5 m i n , one can vary j3 in the definition of 5 to find the maximum value of 6, 
m ax{$min[{T b ) q }}, that is compatible with 5 min : 



8max{Smin[(T b ) q ]} = Maxu arccos 



/r 1 



(3-2) 



The value of 9 max in column (4) is an extreme upper limit because it is the third in a chain 
of bounds. 



• The value of (Tt,) q in table 1 is a lower bound, since NVSS might pick up some extended 
flux that is missed by the peak FIRST measurement and the radio core was actually 
weaker at the time of the NVSS measurement than indicated by the peak flux density, 
i.e., (AF y (mJy)) is underestimated in (2.4). Also, (2.3) is an inequality. From the 
discussion at the end of section 2, the larger (Tf,) 9 , the larger 5 m i n 

• the larger <5 m m, the lower 9 max by (3.2) and the true value of the jet plasma 5 can be 
much larger than 5 min [(T b ) q } 



1 IRAF is the Image Reduction and Analysis Facility, written and supported by the IRAF programming 
group at the National Optical Astronomy Observatories (NO AO) in Tucson, Arizona. 
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• the larger 6 m i n , the lower 9 max by (3.2) and the true value of the jet plasma 5 can be 
much larger than 5 min [(T b ) q ] 

• The actual line of sight to earth satisfies 9 < 9 max by definition. 

Thus the condition for inclusion into table 2, 9 max < 35°, means that the jet is likely to be 
propagating very close to the pole. 



4. The Theory of Polar BALQSOs 



Fig ure 1 is a small modification to Fig. 9 of [Punslyj (1l999bl ). In section 4.5 of iPunsly 



( Il999bl ). it is discussed how a relativistic jet (in red) can coexist nested inside of the bipo- 
lar BAL wind. This jet can emit radio flux that is beamed to within a small angle of the 
polar axis. Within the axisymmetric version of the model, LoBALQSOs exist for lines of 
sight within 15° of the polar axis. The density of the wind is highest nearest the polar 
axis as indicated by the gray shading in figure 1 and LoBALQSOs are viewed closer to the 
polar axis on average than HiBALQSOs. The accretion disk coronal X-rays are screened 
from the BAL wind by the dense base of the jet which provides hydrogen column densities 
of £# ~ 10 25 cm -2 . Similarly, the different lines of sight near the polar axis, through the 
inhomogeneous wind, naturally provide larger Compton scattering columns and more atten- 
uation of the far UV flux from the inner disk and less attenuation of the near UV and optical 
flux from the outer disk, making the spectrum appear red. The more extreme and highly 
polarized LoBALs obviously occ ur in objects wit hout a perfectly axisymmetric distribution 
of dust as discussed in detail in iPunslyl (ll999bl ). Recall that the bipolar BAL wind does 
not preclude the coexistence of a BAL wind from the outer regions of the accretion disk as 
envisioned by iMurray et all (119951 ) . 



Since the LoBALQSOs are viewed closer to the jet axis than the HiBALQSOs in the 
axisymmetric model, geometric arguments imply they should typically have jets with larger 
Doppler factors. Thus, LoBALQSOs should occur in samples of variable BALQSOs at an 
inordinately high rate if the axisymmetric version of the polar model applies to a significant 
subpopulation of the BALQSOs. In our master sample of 116 FIRST, DR5 BALQSOs 
defined in section 3, 69 are HiBALQSOs and 47 are LoBALQSOs. Based on table 1, 4/69 
= 5.80% of HiBALQSOs and 8/47 = 17.02% of the LoBALQSOs have (T b ) q > 5 x 10 12 ° K 
(this condition is equivalent to 9 < 35°). The sample is small, but the fact that in table 1, 
the LoBALQSO likelihood to have a large T b is 2.94 times that of the HiBALQSOs tends 
to support the notion the polar BALQSO model represents a significant subpopulation of 
BALQSOs. 
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Conclusion 



Usin g radio variability arguments, we expanded on the sample of known polar BALQSOs 
begun by IZhou et all (120061 ). In the process, we noted that these radio variable BALQSOs 
have an inordinately large LoBALQSOs subpopulation. It is interesting that these properties 
are expected based on an existing detailed theoretical treatment and modeling of bipolar 
BAL winds JPunslvl E999aB ). It would be informative to continue to monitor the FIRST 



BALQSOs at 1.4 GHz, with matched resolution, in order to find more (T&) 9 > 10 12 K 
BALQSOs and improve the statistical information in table 2. 
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Table 2: Line of Sight to Radio Variable BAL Quasars 



Name 


T b 


^min 


^max 


Type 


(SDSS J) 


(10 12 ° K) 




(degrees) 




081102.91+500724.5 


5.4 


1.7 


34.9 


LoBAL 


081618.99+482328.4 


39.6 


3.4 


17.1 


HiBAL 


081839.00+313100.1 


11.8 


2.3 


26.1 


HiBAL 


082817.25+371853.7 


125.4 


5.0 


11.5 


LoBAL 


093348.37+313335.2 


20.2 


2.7 


21.6 


LoBAL 


104106.05+144417.4 


8.5 


2.0 


29.4 


HiBAL 


113445.83+431858.0 


5.9 


1.9 


33.7 


LoBAL 


134652.72+392411.8 


73.9 


4.2 


13.8 


LoBAL 


142610.59+441124.0 


11.5 


2.3 


26.4 


LoBAL 


145926.33+493136.8 


8.8 


2.1 


29.0 


LoBAL 


155633.77+351757.3 


66.5 


4.0 


14.3 


FeLoBAL 


165543.24+394519.9 


47.5 


3.6 


16.1 


HiBAL 
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Fig. 1. — The bipolar wind model of lPunslvl Jl999al lbl). The gray clouds of the BAL wind 
has an opening angle of ~ 25° - 30°. Nested inside is a relativistic jet (in red). The highest 
BAL wind densities are near the polar axis. Thus, lines of sight near the polar axis, have 
the maximum attenuation of X-ray and ionizing UV radiation from the accretion flow. This 
implies that lines of sight close to the polar axis are more likely to represent LoBALQSOs. 



